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Aim 

To understand how the Fermi-Dirac statistics affect the energy 
of electrons in a metal. 

Objectives 

1. To derive the Fermi-Dirac distribution. 

2. To explain the effect of the exclusion principle on heat 
capacity in a metal. 

3. To explain the meaning of Fermi energy and Fermi 
temperature. 

4. To derive the behaviour of heat capacity at high and low 
temperatures. 
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5.1 Fermi-Dirac statistics 
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Electrons in a metal 



The Drude model was propose in 1900 by Paul Drude to 
explain the transport properties of electrons in a metal. 

This treats the electrons like an ideal gas. 




It explains very well the DC and AC conductivity in metals, the 
Hall effect, and thermal conductivity. 

However, it greatly overestimates the heat capacities. 

http : //en . wikipedia . org/ wiki/Drude_model 
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Why did the electron gas model fail? 



When we derived the Maxwell-Boltzmann distribution for the 
ideal gas, we assumed that it is extremely unlikely for two 
atoms to occupy the same energy level. 
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The reason is that the energy levels are very close together, 
compared to the average energy of the atoms. 

Unfortunately, this is no longer true for electrons in a metal at 
room temperature. 



PHYS393/Hock 



5 



2009-2010 



Recall the equation for the energy of a state: 

h 2 

2 2 2 

Because the electron has a much smaller mass, the energy 
levels are much higher and further apart. 

So the electrons have to occupy all of the lowest levels first, 
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The exclusion principle 



Electrons are not allowed to occupy the same energy states, 



So they have to be stacked up from bottom to top 
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When heated, most of the electrons are stuck - there is no 
space above to move up in energy ! 



Only those few at the very top can. As a result, the heat 
capacities are much smaller than expected of a gas. 
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Fermi-Dirac statistics 



As a result, we cannot use the Maxwell-Boltzmann distribution 
for the ideal gas. 

Fortunately, we can use the same statistical methods that we 
have learnt so far. 

Lets start by looking at the picture of the bundle from another 
angle. 
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We have gi energy states in the bundle. Suppose the bundle 
contains ni electrons. 
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Fermi-Dirac distribution 



According to the exclusion principle, exactly states must be 
filled and exactly gi — ni states unfilled. 

The total number of ways the electrons can be arranged in this 
bundle is therefore: 



For the whole system - i.e. all the bundles - we get 

i 



9i 



% rnKdi ~ rii)\ 

Applying the Lagrange multiplier method again, we would get 

Hi 1 



Qi exp(-Ai - X 2 £i) + 1 
where Ai and A2 are the Lagrange multipliers. 

We have applied the same constraints as before on the particle 
number TV and the energy U. 
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Fermi-Dirac distribution 



Previously, we have used n(s) for the number density, g(e) for 
density of states, and de for the energy interval of the bundle. 

We also know that one of the Langrange multipliers would be 
related to temperature: 

\ 2 = -l/k B T 

The distribution of the electrons may then be written as 

/ x , g(e)de 

niejde — ; r . 

exp(-X 1 +e/k B T) + l 

This is called the Fermi-Dirac distribution function. 
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Fermi-Dirac distribution 



We have obtained the Fermi-Dirac distribution function. 

/ x , g(e)de 

nisjds — ; r . 

exp(-X 1 +s/k B T) + l 

The total number of electrons is fixed and is given by 
integrating the normalisation condition: 

N= f°° g(e)de 

Jo exp(-Ai + e/k B T) + 1 

In principle, we could solve for the Lagrange multiplier Ai. It is 
common practice to express it in terms of \i as follows: 

exp(-Ai) = exp(-n/k B T) 
fi is called the "chemical potential." 
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Occupying the energy states 



For convenience, we define the following function: 



g{e) exp((-^ + e)/k B T) + 1 

Since g(s)ds is the number of states and n(s)ds is the number 
of particles, /(s) would be the fraction of states that are 
occupied. 

So /(s) is called the occupation number. We need to get a feel 
as to what this looks like and how it changes with temperature. 

Lets start with the simplest case: T = OK. If we allow T to 
approach zero, we will find: 

/(s) = l for e < \i 
f (e) = 0 for s > \i 

This means that all states with energy below \i are fully 
occupied. All states with energy above \± are empty. 
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Occupying the energy states 



The graph for f(s) looks like this 



f(e) 



The shape shows that the energy levels are occupied below a 
certain energy, and unoccupied above that. 



This feature is characteristic of the Fermi-Dirac distribution 
that we are studying. 
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The Fermi energy 



At T = OK, The highest energy in the stack of electrons is /x. 
This energy is also called the Fermi energy, E F . 



CD 
i_ 

O 

c 

LU 



We can find the Fermi energy E F by integrating the 
normalisation condition: 



•oo 



N = 2 x n(s)ds 
JO 

The factor of 2 must be added because each energy level can 
be occupied by 2 electrons - spin up and spin down. 



We shall solve this for the Fermi energy E F . 
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The Fermi energy 



In terms of the occupation number, 

So the normalisation condition can be written as: 

roo roo 

N = 2 x / n(s)ds = 2 x / g(s)f(s)ds 

Jo ^ Jo 
We know that at OK, f(s) = 0 for e > /a. So So the integration 
would stop at e = ji\ 

N = 2x [ EF g(e)f(e)de 

Jo 

since \i = E F at OK. 



We also know that at OK, f(s) = 1 for e < /x. So 

rE F 

N = 2 x / g(s)ds 

Jo 
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The Fermi energy 



We have previously derived the density of states: 

g{s) = —^{2msyl 2 

In the topic on ideal gas, this is obtained by counting the 
number of energy states of a particle in a 3-D box. In this 
topic on electrons, we have used the same particle in a box. 
The difference from ideal gas only arises later on, when we 
make different assumptions about the energy levels. 

So the same formula for the density of states can be used for 
both the ideal gas and the electrons. We can therefore 
substitute the formula into the normalisation integral 

rE F 

N = 2 x / g(s)ds 

Jo 

and solve for the Fermi energy. The result is 

H 2 /3tt 2 7V\ 2/3 

E F = 

* 2m \ V 
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5.2 Electronic heat capacity 
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Heat capacity 



We want to use the ideas and formulae that we have developed 
to calculate the heat capacity of electrons. 

To see how to do this, recall that the Drude model would 
predict a heat capacity for that is the same as that of an ideal 
gas: 

C = ^Nk B . 

It is known from experiments that the actual heat capacity of 
the electrons is much smaller. This refers to measurements 
that are done at room temperature. 
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Heat capacity 



We can understand this if we think that the electrons that are 
stacked up to the Fermi energy do not have enough energy to 
jump out of the stack. 



H— 
H- 

This would only be true if the thermal energy is much smaller 
than the Fermi energy. 

In order to find out if this is true, we need to estimate this 
thermal energy. 
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Heat capacity 



We start by assuming that the thermal energy is indeed much 
smaller than the Fermi energy. We shall derive an expression 
for this thermal energy, and then calculate it at room 
temperature to see if the assumption can be justified. 



At a temperature above 0 K, the occupation number 

f(e) = 1 



exp((- M + e)/^T) + l 
would no longer have a sharp step at the Fermi energy. 



f(s) 




If k B T is much smaller than /x, the graph would remain close to 
the step, as if the step has become smoother. 
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Heat capacity 



The "smoothened" slope of the graph tells us that electrons 
just below the Fermi energy = E F ) is excited just above it. 



f(s) 



So we can estimate gain in thermal energy of the excited 
electrons by the width 8 of this slope. 



In order to do this, we take a closer look at the occupation 
number 

f(e) = 1 



exp[(e-/i)/fe B T] + 1 
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Heat capacity 



We are assuming that k B T is much smaller than \±. 



For energy e higher than n by a few times of k B T, The 
exponential function exp[(s — fi)/k B T] would quickly become 
large. 



The occupation number is then approximately 

1 1 

Ke) = 



exp[(s - fi)/k B T] + 1 exp[(s - n)/k B T] 
which is just the exponential function with negative argument 

f(e) = exp[-(e-M)/fe B T]. 
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Heat capacity 



This means that the part of the graph to the right of \i falls off 
exponentially. 




It falls by a fraction of 1/e over an energy range of k B T. 
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Heat capacity 



There is a corresponding behaviour just below the Fermi energy. 

For energy e lower than \± by a few times of k B T, The 
exponential function exp[(s — n)/k B T] would quickly become 
small. 



The occupation number is then approximately 

exp[(s - n)/k B T\ + 1 

where we have used the binomial expansion and kept only the 
first order term. 
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Heat capacity 



This means that the part of the graph to the left of /x tends to 
the line f(e) = 1 exponentially. 



It reaches within f(e) = 1 by a fraction of 1/e, over an energy 
range of k B T. 

This means that it is the electrons within this energy range 
that is excited. So the thermal energy of the excited electrons 
is of the order of k B T. 
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f(s) 




Heat capacity 



At this point, we should justify our assumption that k B T at 
room temerature is much smaller than the Fermi energy 
IJL = e Fi which is defined at 0 K. 

We shall take sodium metal as an example, and calculate kgT 
and Ep for this metal. 

In sodium, each atom has one valence electron. This electron 
is mobile and forms the electron gas that we are talking about. 
In order to calculate the Fermi energy, we need the number 
density N/V. We can calculate this from the following data: 

density = 0.97 g cm -3 

relative atomic mass = 23.0 

So the volume for one mole of atoms is 

23^-0.97 = 23.71 cm 3 , 
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Heat capacity 



and the number density is 

TV 

— = Na 23.71 

V 

where Nj± is Avogadro's number. This gives an answer of 
2.54 x 10 28 m~ 3 . 

Using the Fermi energy formula, 

h 2 f3n 2 N\ 2/3 
E F = 

where m is the mass of the electron, we can find that the 
Fermi energy is 3.16 eV. 
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Heat capacity 



In contrast, at room temperature 298 K, we can calculate that 

1 

k B T « — eV. 
B 40 

This is about 120 times smaller than the Fermi energy. 

We can repeat this for other typical metals, and we would get 
similar answers. 

This justifies our assumption that at room temperature k B T is 
much smaller than the Fermi energy. 

We are now a step closer to estimating the electronic heat 
capacity. Next, we need to understand the behaviour of the 
excited electrons. 
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Heat capacity 



When temperature increases above 0 K, the step in the 
Fermi-Dirac distribution becomes smoother as electrons just 
below the Fermi energy are excited above it. 




Notice that the "tail end" of the distribution - to the right - 
looks exponential. This is because there are relatively few 
electrons above the Fermi energy. So these can behave like the 
ideal gas and approximately obey the Maxwell-Boltzmann 
distribution. 
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Heat capacity 



It is clear from the graph that for a small increase in 
temperature, only electrons close to the Fermi energy are 
excited. Most of the electrons are below the Fermi energy and 
are not excited at all. Since they cannot be excited, these 
electrons would not contribute to the heat capacity . 




So it is mainly the electrons close to the Fermi energy that 
would contribute to the heat capacity. So we can use these to 
estimate the heat capacity and ignore the rest. 
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Heat capacity 



As we have seen, these electrons close to the Fermi energy 
behave like the ideal gas. We know that the energy of an ideal 
gas is 

where N\ is the number of particles in the ideal gas. 

In the case of the electrons, N\ should refer to the number of 
electrons above the Fermi energy, and not the total number of 
electrons. We can estimate this number as follows. 
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Heat capacity 



The electrons that do get excited are in the small energy range 
of order kgT from the Fermi energy E F . 




1-exp[(s-^L)/k B T] 



So we can estimate N± with the number of electrons that are 
within the energy interval of 



de = k B T 



from the Fermi energy. 
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Heat capacity 



f(s) 




1-exp[(s-jLi)/k B T] 



We know that the number of particles in a given energy 
interval is 

n(s)ds = 2g(s)f(s)ds 

where the factor of 2 again comes from the spin states of the 
electrons. 



At OK, the energy states below E F are fully occupied, i.e. 
f( £ ) = i. So the number would given by 

TVi 2g(E F )k B T 
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Heat capacity 



From the definitions of the density of state and the Fermi 
energy, it can be shown (see exercises) that 

g(E F ) = 



4E F 

where TV is the total number of electrons. This gives 

/ 37V \ 3 k B T 

Substituting into the energy for ideal gas 

Ui = \N x k B T 
we get the energy for the excited electrons 



2 V 2 Ep J 4 E F 

Differentiating with respect to T, we get the electronic heat 
capacity 

C = 9 -Nk B k -^ 

2 _ E F 
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Heat capacity 



We summarise what we have learnt about the heat capacity: 

Some states below Fermi level becomes empty, and some 
states above becomes occupied. 



k B T 



For a temperature T, an electron that gets excited would be 
able to gain on average an energy of about 3k B T/2. 
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Heat capacity 



1. Only electrons within k B T of the Fermi level can jump up. 
Those further down cannot jump because the states above are 
still mostly occupied. 

2. The number of electrons that can get excited is 3 2E^ N ■ 
This is a fraction 3 Jf^ of the total. 

3. Each would gain an energy of about 3/c#T/2. So the 
increase in total energy is obtained by multiplying the number 
of excited electrons by this energy 

AU « 31 ^N x 3k B T/2 = -Nks 1 ^. 
2E F 1 4 E F 

4. Differentiating with respect to T, we get the electronic heat 
capacity 

C = 9 -Nk B k * T 



2 ~ E 



F 
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Heat capacity of a metal 



We have obtained the electronic heat capacity 

More detailed calculations show that the factor of 9/2 should 
really be 7r 2 /2: 

C = ^Nk B k * T 



2 E F 

This equation is often written in the form 

C = —Nk B 



2 T F 

where T F = E F /k B is called the Fermi temperature. 



PHYS393/Hock 



37 



2009-2010 



Heat capacity of a metal 

Notice that the heat capacity is directly proportional to T. 

tt 2 T 

C = —Nk B — 

This is often written in the form 

C = tT. 

There is another contribution to the heat capacity. This comes 
from the vibrations of the atoms, and it is proportional to T 3 . 

We could imagine writing the total heat capacity in the form: 

c v = 1 T + AT 3 . 

We can measure this heat capacity to check if the formula is 
correct. Suppose that we have obtained a table of values for T 
and c v . To check if the formula is correct, we can rewrite it in 
this form: 

c v /T = 1 + AT 2 
If we plot cy/T against T 2 , we should get a straight line. 
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Heat capacity of a metal 



In 1955, William Corak and his fellow co-workers measured the 
heat capacities of copper, gold and silver from IK to 5K in 
their laboratory in Pittsburgh ... 




FlG. 6. Atomic heat of silver. 



W. S. Corak, et al, Physical Review, vol. 98 (1955) pp. 1699-1707 

and got the straight lines. This shows that the predictions of 
the Fermi-Dirac statistics are correct. 
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5.3 Effective mass, and a few other things 
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Heat capacity of a metal 



From the formula, 

c v /T = 1 + AT 2 

we know that we can find 7 from the y-intercept of the graph. 
The following table shows the values for a few metals. 



FREE ELECTRON 7 MEASURED y RATIO 0 

ELEMENT ^ 1Q _ 4 ^^^-1^-2) (m*fm) 

Li 1.8 4,2 2.3 

Na 2.6 3.5 1.3 

K 4.0 4.7 1.2 



The second column contains the values of 7 from the formula: 

C = *-Nk B k -^ = 7 T 
2 B E F 1 

The third column contains the values actually measured. They 
are obviously different. 
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Heat capacity of a metal 



Does this mean that our theory is wrong? 

Yes and no. The theory tells us that the electronic heat 
capacity is proportional to T. The measurement shows that 
there is indeed such a contribution. 

So we are not completely wrong. Perhaps the theory needs 
refining. We can be optimistic and go back and try and 
understand what we have missed. 

Recall that we have started with a particle in a 3-D box and 
calculated the energy levels. Then we just fill these up with 
electrons and calculated the heat capacity. All we have is a gas 
of electrons in empty space. 
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Heat capacity of a metal 



But what about the atoms? The "empty space" is really filled 
with atoms. The electrons must surely interact with the atoms. 
This then is the reason for the difference between theory and 
measured 7. 

However, the measured C is proportional to T. This agrees 
with theory, and it should mean something. One possibility is 
that, for some reason, the electron interacts only weakly with 
the atoms. This idea has been shown to be correct by other 
types of measurements. 

According to this idea, the behaviour of the electrons in the 
presence of the atoms is essentially the same. The difference is 
that the interaction with the atoms make the electrons behave 
as if they have a different mass. 
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Heat capacity of a metal 



If we imagine that the electrons in a metal has a different 
effective mass ra* than its natural mass, we can "explain" the 
difference in 7. 

According to the formula for heat capacity, calculated 7 is 
given by 



which is in turn inversely proportional to the mass m of the 
electron. 

So 7 is directly proportional to the mass. 




This is inversely proportional to the Fermi energy 
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Heat capacity of a metal 



This means that we can get the effective mass if we divide the 
measured 7 (in the third column) by the calculated 7 (in the 
second column). 

This is shown in the last column of the table. 



FREE ELECTRON 7 MEASURED f RATIO 0 

ELEMENT ^ 1Q _ 4 ^^^-1^-2) (m*fm) 

Li 1.8 4,2 23 

Na 2.6 3.5 13 

K 4.0 4.7 1.2 
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Behaviour at high temperatures 



At very high temperatures, all the electrons could get excited. 
Then they would start to behave like an ideal gas. There are 
many more energy state they can reach and they are much less 
likely to be forced towards the same energy states. 

This would happen only if the electrons are excited far above 
the Fermi energy. Therefore we can use Fermi temperature as 
a reference. A temperature would be high if it is high compared 
to the Fermi temperature. 

At the high temperature, the heat capacity would therefore 
change to that of the ideal gas. So instead of 




T 



which is very small, it would become 



which is much larger 
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A note on the chemical potential 



We have seen many times that the chemical potential jjl is 
equated to the Fermi energy E F . It is important to note that 
this is true only at 0 K. 

Recall the chemical potential n is determined by the number of 
particles: 



which is why \i is the Fermi energy then. 

At higher temperature, \i would change. But for temperatures 
well below the Fermi temperature, as with electrons in metals 
at room temperature, we may assume that \i is still quite close 
to the Fermi energy. 
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At T = 0 K, the integral is just 




A note on fermions 



Note that electrons belong to a larger family of particles called 
fermions. 

A fermion is a particle with a half integer spin. Other examples 
are 

the proton, 
the neutron, 

the helium-3 ( 3 He) atom and 
the oxygen-13 ( 13 0) atom. 

All fermions are known to show the kind of properties we have 
seen for electrons. They obey the Fermi-Dirac statistics. 
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What we have learnt so far 



1. Electrons follow the Fermi-Dirac distribution: 

n(s)ds = . 

exp[(e-fi)/k B T] + l 

2. At 0 K, the electrons fill up the energy states from the 
lowest level. The highest energy of the electrons is called the 
Fermi energy: 

h 2 /3tt 2 7V\ 2/3 



Ep — 



2m \ V 



3. When thermal energy k B T is small compared to the Fermi 
energy, only electrons near the Fermi energy would get excited 
So only these would contribute to the electronic heat capacity: 

tt 2 T 
C = —Nk B 



2 T F 

4. The Fermi temperature is defined as 

k B 
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5.4 Exercises 
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Exercises 



Exercise 1 



Show that the density of states at the Fermi energy E F is given 
by 

9{E F ) = 



4E F 

where TV is the number of electrons. 
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Exercises 



We can make use of the following relations for the density of 
states. 

Relation with the total number of states: 

, x G{e) 

g(£) = — ■ 

as 

Formula for the total number of states: 

G(e) = ^(2 W£ )^ 

SO 

G(e) oc s 3 / 2 . 
Relation with the Fermi energy and TV: 

rEp 

N = 2 g(e)de. 

Jo 
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Exercises 



From 



and 



we have 



From 



we have 



, x G(s) 
as 



N = 2 I g(e)de. 

Jo 



N = 2G(E F ) 



G(e) <x e 3 / 2 . 



sOO 3 



GO) 2e 
At the Fermi energy Ep, we get ... 
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Exercises 



_ 3G(E F ) 
Substituting the above result for TV: 



we get 



N = 2G(E F ), 



9(E F ) = 



AE 



F 
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Exercises 



Exercise 2 

The number of ways of arranging fermions is given by 



i n iK9i ~ rii)\ 

Use this formula, along with the conservation of particle 
number and total energy, to derive the Fermi-Dirac distribution 
function: 

n% 1 

Qi exp(-Ai - X 2 €i) + 1 

where Ai and X 2 are the Lagrange multipliers. 
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Exercises 



We need to maximise the Lagrange function In Q + XiN + X 2 U. 

Using Stirling's formula 

In TV! « 7Vln N -N 

we have 

Differentiating with respect to r^, 

— = - In rii + ln(& - = In ^ * 

The total number is given by 

N = J2 n 

i 

and the total energy by 

U = n,.s 



Differentiating these ... 
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Exercises 



We get 

ON 



= 1 

dni 

and 



dU 
dni 



Then differentiating the Lagrange function, we get 

ln (^r) + Al + A2£! = 0 

Rearranging, we get the Fermi-Dirac distribution function 

Hi 1 



Qi exp(-Ai - X 2 si) + 1 
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Exercises 



Exercise 3 

Show that the mean energy of electrons at absolute zero, e is 
3Ep/5, where E F is the Fermi energy. 
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Exercises 



The total energy is 

POO 

U = n(s)eds. 
JO 

For the Fermi-Dirac distribution, 



roo 

U = 2 g{e)f{e)ede. 

Jo 



At 0 K, this would be 



rE F 

U = 2 g(s)sds. 

Jo 



Recall the density of states 



, 4-m7rV . .1 /o 
g(e) = -^(2ms) 1 / 2 



Substituting into the integral, 



f E F A-m7TV , N q/o 

Jo h 3 y J 
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Exercises 



Integrating, 

U = 2 



4m7rV (2m)^ 2 E 5 J 2 



h3 v-» v 

To relate this to N, 

rE F 

N = 2 g{e)de. 

Jo 

we substitute the density of states 

Jo h 3 y J 



and integrate: 

TV = 2 



Wy (2m) 3 / 22 ^ 2 



^ 3 ^.. VJ ^ F 

The average energy is the energy per particle 

- U 3^ 
e = — = -E F 
N 5 * 
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Exercises 



Exercise 4 

Metallic hydrogen can be produced by compressing liquid 
hydrogen to extremely high pressure. This transition has been 
observed when the density reaches 0.32 x 10 6 mole m -3 , at a 
temperature of 3000 K. It is thought that in the metallic phase 
the liquid comprises ions and electrons. Were the electrons 
in this liquid metal obeying classical or Fermi-Dirac statistics? 
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We need to compare the 3000 K temperature to the Fermi 
temperature T F . 

If the temperature is much higher than T F , the electrons would 
obey classical statistics and behave as an ideal gas. If the 
temperature is much smaller than T F , the electrons would obey 
Fermi-Dirac statistics. 

The Fermi temperature is related to the Fermi energy by 



and the Fermi energy is given by 




2/3 



where m is the mass of the electron. 
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The hydrogen density 0.32 x 10° mole m - ^ is also the electron 
density, since each molecule H2 converted to gives one 
electron. The number density is obtained by multiplying 
Avogadro's constant: 

— = (0.32 x 10 6 ) x (6.022 x 10 23 ). 

Substituting this and the other constants into the above 
formulae, we get the Fermi temperature 

T F = 141,000*:. 

The temperature of 3000 K is therefore much lower than the 
Fermi temperature. So the electrons in metallic hydrogen 
would obey Fermi-Dirac statistics, in the same way as normal 
metals at room temperature. 
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Exercise 5 

One possible means of obtaining fusion energy is to implode 
spherical capsules containing heavy hydrogen by irradiating 
them with high power lasers. For fusion to occur the implosion 
core needs to reach electron temperatures of the order of ~ 1 
eV, and the electron density needs to be of the order of 
10 33 m -3 . Do we need to use Fermi-Dirac statistics to describe 
the electrons, or are Maxwell-Boltzmann statistics sufficient? 
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We need to compare the 1 eV electron temperature to the 
Fermi energy E F . The Fermi energy is given by 



The number density N/V is lCr d m J . Substituting into the 
above formula, we get 

E F = 3650 eV. 

The electron temperature of 1 eV is much smaller than this. 

Therefore the electrons would be described by Fermi-Dirac 
statistics. They would behave like electrons in metal at room 
temperature, and not as an ideal gas. 




2/3 
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Note 

Notice that the terms energy and temperature have been used 
interchangeably. 

A temperature of E = 1 eV, means that the temperature 
would be 

T = E/k B 

in Kelvin. 
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